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Abstract 

Pseudonaja textilis is a widespread and common snake in eastern parts of Australia, but its distribution in New Guinea is 
poorly understood, and the origin of the New Guinea populations and its timing have been the subject of much specula¬ 
tion. Phylogenetic analysis of mitochondrial DNA sequences from three New Guinea populations ofP textilis indicates 
that New Guinea was colonised from two independent eastern and western migration routes most likely in the Pleis¬ 
tocene. One dispersal event from northern Queensland led to the populations in eastern New Guinea (Milne Bay, Oro and 
Central Provinces, Papua New Guinea), whereas another, from Arnhem Land to central southern New Guinea, led to the 
populations from the Merauke area, Indonesian Papua. The results are consistent with the effects of Pleistocene sea level 
changes on the physical geography of Australasia, and are thus suggestive of a natural rather than anthropogenic origin 
of the New Guinea populations. The taxonomic status of the New Guinean populations is discussed. 

Key words: Pseudonaja textilis , New Guinea, Australia, phylogeography, Pleistocene, sea level changes, Lake Carpen¬ 
taria, taxonomy, mitochondrial DNA 


Introduction 

In recent years, the advent of molecular markers, and particularly the use of mitochondrial DNA-based phylo- 
geographical approaches, has revolutionised our understanding of the biogeographical history of animal spe¬ 
cies (Avise, 2000). One region that has received increasing attention in this regard is the Australasian region, 
including the faunal relationships between Australia and New Guinea. While earlier studies focussed prima¬ 
rily on mammals (e.g., Aplin et al. , 1993; Pacey et ai, 2001), reptiles, including snakes, have been studied 
extensively in recent years (Harvey et al. 2000; Rawlings & Donnellan, 2003; Rawlings et al., 2004; Kuch et 
al., 2005; Wiister et al., 2005). The results have revealed a variety of biogeographical patterns, including older 
exchanges resulting in endemic New Guinean species that are clearly distinct from their Australian sister taxa 
(e.g., Acanthophis laevis, Pseudechispapuanus and P. rossignolii —Kuch et al., 2005; Wiister et al., 2005), to 
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lower levels of divergence, where Australian and New Guinean sister groups are nested within a wider 
intraspecific phylogeny (e.g., Acanthophis rugosus —Wiister et al., 2005; Morelia viridis —Rawlings & Don- 
nellan, 2003; Morelia amethistina complex—Harvey et al., 2000), and to zero or minimal differentiation (e.g., 
Oxyuranus scutellatus —Wiister et al., 2005; Doughty et al., 2007) 

Among the elapid snakes found on both sides of the Torres Strait, the origins of the New Guinea popula¬ 
tions of the genus Pseudonaja have remained the most contentious. The genus Pseudonaja is widespread in 
Australia, where it is represented by seven species (Cogger, 2000), although several analyses have revealed 
the likely presence of cryptic species within what is currently referred to as Pseudonaja nuchalis Gunther 
(Mengden, 1985; Skinner et al., 2005). A single species, P. textilis, inhabits both Australia (where it is wide¬ 
spread and common in eastern parts, but more patchy and sporadic in the centre) and parts of the island of 
New Guinea. 

The earliest records of P. textilis in PNG were a female from Menapi, Cape Vogel, and a male from 
Baiawa, Moi Biri Bay (part of Collingwood Bay), northern Milne Bay District, collected during the 4 th Arch¬ 
bold Expedition to New Guinea in 1953 (AMNH 73949 & 73959 respectively) (McDowell, 1967). McDowell 
further remarked that the two specimens resemble Queensland and North Australian Pseudonaja in all 
respects except one, the presence of 12, as opposed to 9-11, solid maxillary teeth. 

Eric Worrell (1961) received a specimen at the Australian Reptile Park, Gosford, shipped with Pseudechis 
papuanus from locations in ‘Papua’, which then referred to the southern part of present-day Papua New 
Guinea (Campbell, 1969). However, in 1961 and 1963 he recorded the species’ New Guinea distribution erro¬ 
neously as “Highlands” (Worrell, 1961). No records or specimens remain. 

Slater (1968) included P. textilis in the second edition of his small handbook, referring to specimens being 
collected around Dogura, Milne Bay District. His source was probably the Dogura specimen which arrived in 
Port Moresby in 1965 and which was subsequently identified as P textilis by Cogger at the Australian 
Museum (Campbell, 1969). It is likely that Slater was unaware of the two Archbold specimens. Campbell 
(1969:26) mentioned the Dogura specimen and also reported a further two specimens from Dogura (AM 
R25774-5), and the head of a specimen from Embogo Hospital (AM R25880), the first record from Oro (then 
Northern District), collected around the same time and identified as P. textilis by Cogger. The National 
Museum and Art Gallery of Papua New Guinea also contains a specimen from Dogura (as Dogua) (NMPNG 
22295—old ace. no. 10419). 

The records from the northern coast of eastern PNG, in Milne Bay and Oro Provinces, remained the only 
confirmed records until the 1990s. This apparently atypical distribution for a savanna-dwelling elapid with 
Australian affinities (the remainder of which are restricted to the savannas fringing the southern side of New 
Guinea —Oxyuranus scutellatus, Pseudechis papuanus, P. rossignolii, Demansia vestigiata ), and the lack of 
records of the species in apparently suitable habitat in Papua New Guinea’s Central and Western Provinces, 
led Slater (1968) to suggest that the species’ presence in New Guinea may have been the result of an acciden¬ 
tal introduction either in military equipment imported during or after World War II, or in post-war agricultural 
supplies, a hypothesis commented on by subsequent workers (e.g., O’Shea, 1990, 1996; Lalloo, 1994; Kuch & 
Yuwono, 2002; Williams & Wiister, 2005). 

The presence of Pseudonaja in the Central Province savannas has long been suspected. Early clues to its 
presence came from five positive ELISA results from snakebite victims (Lalloo, 1994; O’Shea, 1996) from a 
number of scattered localities in Central Province and the National Capital District (NCD). The first specimen 
from the southern savannas, NMPNG 24583 (old ace. no. 18937) from 6-Mile (NCD) was discovered by the 
authors in 1996, in the NMPNG collection, mislabelled as Oxyuranus scutellatus. The head of a large speci¬ 
men from Hisiu (Hiri Sub-Province, Central Province) is also represented (NMPNG 25305). Four additional 
cases of envenoming by P. textilis were detected using immunological assays among patients admitted to Port 
Moresby General Hospital from Central Province during 2006 (D.J. Williams, in prep.). In combination, these 
scattered specimens, immunoassay results and snakebite records suggest that P. textilis is broadly distributed 
throughout the southern savannas of Central Province. 
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Pseudonaja textilis remains a rarely collected snake in PNG, but the same cannot be said of the population 
in West Papua. Currie (2000) believed that Pseudonaja sp. could inhabit the region but was unaware that J.S. 
Keogh had placed two specimens from Merauke in the Australian Museum in 1995 (AM R147652, R147659). 
Kuch and Yuwono (2002) reported that since 1993 numerous specimens have been captured in the vicinity of 
Merauke. Pseudonaja textilis could be expected to inhabit the neighbouring savannas of Western Province, 
PNG, but to date no positive record or specimen exists. Based on the recent records suggesting a widespread 
distribution of the species in New Guinea and the abundance ofR textilis in the vicinity of Merauke, Kuch and 
Yuwono (2002) concluded that the hypothesis of an origin through anthropogenic introduction was an 
unlikely explanation for the presence of P. textilis on the island. 

The taxonomic status of the New Guinea populations of Pseudonaja has remained unclear. McDowell 
(1967) assigned his material to Pseudonaja textilis, but other authors (e.g., O’Shea, 1996; Kuch & Yuwono, 
2002) considered the status of the form as unresolved, referring to the New Guinea populations as “ Pseud¬ 
onaja cf. textilis ”. As in the case of many other Australasian snakes, the nomenclature of the genus Pseud¬ 
onaja has been confused by largely evidence-free contributions from a small minority of amateur 
herpetologists (Williams et al., 2006b; see Skinner et al., 2005, for a discussion of the genus-level classifica¬ 
tion of Pseudonaja). Using the Archbold specimens described by McDowell (1967) as type material, Hoser 
(2003) described the New Guinea populations of P. textilis as a separate subspecies, P. textilis pughi, empha¬ 
sising that this would include all New Guinea populations of Pseudonaja, including those from Merauke, 
West Papua. His diagnosis was based primarily on the dentitional differences between Milne Bay and northern 
Queensland specimens reported by McDowell (1967), as well as anecdotal evidence of slight differences in 
colour. 

The recent study of the phylogeny and phylogeography of Pseudonaja by Skinner et al. (2005) included a 
large number of geographically widespread samples of P. textilis from Australia, and Keogh’s two specimens 
from Merauke (Indonesian West Papua), and revealed considerable phylogeographic structuring within P. tex¬ 
tilis. The specimens from Merauke were found to nest among a set of haplotypes from northern South Austra¬ 
lia, western Queensland and the Northern Territory. However, representatives of the apparently allopatric 
populations from eastern parts of New Guinea were not included in that study. 

Understanding the origin and taxonomic status of the New Guinea populations of P. textilis is especially 
important from the point of view of snakebite management in New Guinea (Williams & Wiister, 2005). In 
Australia, Pseudonaja spp. are the elapid snakes that have adapted best to anthropogenic habitat changes, and 
as a result, they have become the most common cause of snakebites and snakebite fatalities (Sutherland, 1992; 
White, 1998) in Australia. At present, bites by P. textilis are an apparently uncommon cause of snakebite mor¬ 
bidity and mortality in New Guinea (Lalloo et al., 1994; D.J. Williams, in prep.). Understanding the origin 
and taxonomic status of the New Guinea populations of P. textilis, and monitoring their abundance and distri¬ 
bution, is essential to predict likely future trends in the medical importance of this species, and thus antivenom 
needs for New Guinea, especially as human impacts on vegetation and habitats may create conditions favour¬ 
ing this species over other native elapids. In particular, it is essential to determine whether the late discovery 
of the species and the currently accumulating records of its presence reflect rarity in a native species restricted 
to a relict distribution that has been neglected by collectors, or whether this is a reflection of a recent introduc¬ 
tion and relatively rapid range expansion. The latter case predicts a rapid future increase in the medical impor¬ 
tance of this species in New Guinea, whereas in the former case, studies on the effectiveness of existing 
antivenoms against Pseudonaja for the New Guinea populations may be required (Fry et al., 2003). 

During 2006, three Pseudonaja were collected by CCA from Wamawamana and environs on the north 
coast of Milne Bay Province, and one of us (MOS) observed six additional specimens and captured three of 
them in an oil-palm plantation at Heropa, between Doboduru and Embogo on the Buna road, Oro Province, 
Papua New Guinea. These recent new specimens, and the combination of the mitochondrial DNA sequences 
obtained from these specimens with the published data of Skinner et al. (2005), provide an opportunity to 
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investigate the origin and status of these populations, with the aim of resolving several questions, in particu¬ 
lar: 

• Do all New Guinea P. textilis descend from a single colonization event? 

• Is the phylogenetic nesting of the New Guinea haplotypes consistent with natural colonisation from Aus¬ 
tralia, or does it indicate human intervention, or possibly a combination of both? 

• Is there evidence that the New Guinea populations of P. textilis represent one or more distinct taxa from 
those in Australia? 


Material and methods 

Blood or tissue samples were obtained from five PNG specimens (AVRU-UPNG live collection PT001- 
PT003 and LSUMZ 90636-37), as well as from an additional specimen from Merauke, Papua, Indonesia (gift 
from Ulrich Kuch). Total DNA was extracted by standard methods (Sambrook et al., 1989). 

Laboratory work was conducted in two independent labs. A ~ 900 base pair fragment of the mitochondrial 
NADH dehydrogenase subunit 4 (ND4) gene and the adjoining tRNA-His and tRNA-Leu genes was ampli¬ 
fied using either the primers ND4 (Arevalo et al., 1994) and HIS12763v (5’-TTC TAT CAC TTG GAT TTG 
CAC CA-3’, Sylvain Ursenbacher, pers. comm.) or the primers used by Skinner et al. (2005) (ND4 and 
M246). The two PCR protocols were used in each lab: (1) 20 pi reactions were carried out using 18pl of 2x 
ReddyMix™ PCR Mastermix (Abgene™, catalogue no. AB-0575-DC-LD), consisting of 1.25 units of Ther¬ 
moprime Plus DNA polymerase, 75mM Tris-HCL (pH 8.8 @ 25°C), 20mM (NH 4 ) 2 S0 4 , 1.5 mM of MgCl 2 , 
and 0.01% (v/v) Tween® 20, 0.2 mM of each dNTP and precipitant red dye for electrophoresis. (2) 25 pi reac¬ 
tions were carried out using 0.2 mmol each dNTP, 5 pmol each primer, 0.5 units of Taq polymerase (New 
England BioLabs Inc., Ipswich, MA), 2.5 pL of 10X Buffer and 18.9 pL nuclease-free, ultrapure H 2 0 in 25 
pL. 

The amplification protocol comprised denaturation at 94°C for 2 minutes, followed by 35 cycles of dena¬ 
turing for 45 seconds at 92°C, followed by 1 minute of annealing at 50-55°C, one and a half minutes exten¬ 
sion at 72° C. The reaction ended with one longer extension phase at 72° C for 5 minutes. Sequencing was 
carried out either by Macrogen (Seoul, S. Korea—http://dna.macrogen.com), using the ND4 primer or PCR 
amplicons were purified with ExoSAP-IT® (USB Corp., Cleveland, OH), cycle sequenced for both ND4 and 
M246 primers using BigDye® v3.1 and run on an ABI-3100 Genetic Analyzer (Applied Biosystems, Foster 
City, CA). 

Newly generated sequences were aligned with those of Skinner et al. (2005), GenBank accession numbers 
DQ098598-DQ098645. To test for the presence of nuclear pseudogenes (Zhang & Hewitt, 1996), we trans¬ 
lated the DNA sequences into amino acid sequences in MEGA 2.1 (Kumar et al., 2001) to check for prema¬ 
ture stop or nonsense codons or frameshifts. 

Before phylogenetic analysis, only single copies of each haplotype were retained. We noted no differences 
between the GenBank sequences of haplotypes 78 and 85, and haplotypes 96, 99 and 101 from Skinner et al. 
(2005), and therefore did not include them separately. We tested for the presence of a significant phylogenetic 
signal by examining the skewness of the distribution of the lengths of 10 6 random trees (Hillis & Huelsenbeck, 
1992). For phylogenetic analysis, we used maximum parsimony (MP) and Bayesian inference (BI) methods. 
MP analysis was carried out using the software PAUP* 4.0bl0 (Swofford, 2002). For MP analysis, we per¬ 
formed an unweighted heuristic search, using TBR swapping and 10,000 random addition sequence repli¬ 
cates. Internal support for different nodes was estimated using non-parametric bootstrap searching 
(Felsenstein, 1985), using 1000 bootstrap replicates and heuristic searches with 10 random addition sequence 
replicates each. We also estimated Bremer support (Bremer, 1994) by repeating the heuristic searches while 
retaining successively longer trees, and determining how many steps were necessary to break up each clade. 
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Trees were rooted using outgroup sequences from Pseudonaja affinis Gunther (DQ098468), P inframacula 
(Waite) (DQ098467) and P ingrami (Boulenger) (DQ098488). 

For phylogenetic analysis using BI, we used MrBayes 3.1 (Ronquist & Fluelsenbeck, 2003). There is 
increasing evidence that complex models of sequence evolution can extract additional phylogenetic signal 
from data, especially where saturation of base pair substitutions is commonplace (Castoe et a!., 2004). There¬ 
fore, we used different models of sequence evolution for biologically relevant partitions of our data. In the 
case of protein coding mitochondrial genes, the most relevant partitions are first, second and third codon posi¬ 
tions, which are known to display different patterns of sequence evolution. We therefore partitioned our data 
into four separate data partitions, namely first, second and third codon positions separately for ND4, and the 
tRNA sequences as a single partition. To identify the most appropriate models of sequence evolution for each 
data partition, we used MrModeltest 2.2 (Nylander, 2004), and selected the model favoured under the Akaike 
Information Criterion for each category in our Bayesian analysis. Since MrBayes can only use a single out¬ 
group taxon, P. ingrami was specified as the sole outgroup for all BI analyses. We ran the analysis for 4 x 10 6 
generations using 4 simultaneous independent runs initiated with different random starting trees, and sam¬ 
pling every 200 generations. Plots of InL against generation were inspected to determine the bum-in period, 
and trees generated prior to the completion of bum-in were discarded. As a safety margin, we discarded the 
first 10 6 generations. 

To test whether the data reject the monophyly of the haplotypes found in either the New Guinean or the 
Australian populations of P. textilis with statistical significance, we used PAUP* to build the most parsimoni¬ 
ous trees constrained to be consistent with those hypothesis. We then used Wilcoxon signed-ranks tests (Tem¬ 
pleton, 1983, implemented in PAUP*) to test for the significance of differences in tree length between the 
most parsimonious trees and the constrained tree. We used the Shimodaira-Hasegawa test (SH test—Shimo- 
daira & Hasegawa, 1999) to perform the analogous analysis under the maximum likelihood (ML) criterion. To 
do this, we carried out a heuristic ML search in PAUP*, using TBR branch swapping and 10 random addition 
sequence replicates, and a model of sequence evolution obtained for the unpartitioned data with MrModeltest. 
The ML search was then repeated while constraining the New Guinean or the Australian haplotypes to be 
monophyletic. The constraint trees were compared to the optimal ML tree using the RELL option of the SH 
test in PAUP*. Finally, in the case of Bayesian inference analysis, we used PAUP* to filter all trees obtained 
after completion of the bum-in phase to retain only those consistent with the alternative constraint trees, and 
considered the alternative hypothesis as rejected if it was supported by less than 5% of Bayesian trees. 


Results 

We obtained sequences corresponding to those of Skinner et al. (2005) from the new material from New 
Guinea (GenBank accession numbers for the Milne Bay and Oro specimens: EU258951-EU258952; our addi¬ 
tional Merauke sample yielded an identical sequence to Skinner et aids Merauke haplotype). There were no 
indels, unexpected stop codons or frameshifts in the protein coding parts of our sequences, leading us to con¬ 
clude that they represent genuine mitochondrial DNA (Zhang & Hewitt, 1996). Of 767 included base pairs, 
148 were variable and 74 parsimony-informative. The 10 6 random trees generated a skewness statistic of gl= 
-0.662122, rejecting the null hypothesis that the data contain no significant phylogenetic signal (P < 0.001— 
Hillis & Huelsenbeck, 1992). 

Parsimony analysis resulted in 20 equally most parsimonious trees of 202 steps (c.i. = 0.7673, r.i. = 
0.823). MP and BI analysis revealed essentially similar phylogenetic trees (Fig. 1). The results are largely 
identical to those of Skinner et al. (2005). Both analyses reveal three distinct clades of haplotypes within P. 
textilis (Fig. 2): a robustly supported southeastern clade 1, with samples from southeastern South Australia, 
Victoria, and southern and eastern New South Wales, a weakly supported northeastern Australian clade 2, 


ORIGIN OF PSEUDONAJA IN NEW GUINEA 


Zootaxa 1703 © 2008 Magnolia Press ■ 51 



with samples from eastern South Australia, central and northern New South Wales and eastern Queensland, 
and a robustly supported central Australian clade 3, with samples from scattered localities in the Northern Ter¬ 
ritory, northern South Australia and western Queensland. The main difference between our results and those 
of Skinner et al. (2005) is the placement of the haplotype from Malanda, northern Queensland, which clus¬ 
tered with clade 3 in Skinner’s Bayesian analysis (with low support), but with clade 2 (and stronger support) 
in his parsimony tree. Our data consistently place it at the base of clade 2, albeit with low support values. The 
samples from New Guinea are diphyletic, nesting separately within two of the three clades: the samples from 
Milne Bay and Oro Provinces, eastern Papua New Guinea, nest within clade 2, whereas the Merauke haplo¬ 
type nests within clade 3. 



FIGURE 1 . Bayesian inference tree for the included haplotypes of Pseudonaja textilis. New Guinea haplotypes are in 
bold type and enlarged font, all other samples are Australian. Branch support values are Bayesian posterior probabilities, 
parsimony bootstrap and Bremer support. Only values of over 0.5 (Bayesian) and 50% bootstrap support (MP) are given, 
dashes indicate support below that level. Taxon label information includes GenBank accession number, locality, state/ 
country (NSW = New South Wales, SA = South Australia, Qld = Queensland, NT = Northern Territory, PNG = Papua 
New Guinea), and, where applicable, the haplotype number from Skinner et al. (2005). 
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FIGURE 2. Distribution of sampling sites, haplotype clades and hypothesised colonisation routes of Pseudonaja textilis 
into New Guinea. + indicates populations from Central Province, PNG, that were not sampled in this study, but are pre¬ 
dicted to nest with the eastern New Guinea populations. Light grey shading designates likely areas of dry land at times 
with sea levels more than 75 metres below current levels (modified from Voris, 2000). This is only shown for the Sahul 
Shelf, not for other parts of the Indonesian islands or Australia. Dark grey shading indicates the approximate location of 
Lake Carpentaria (Keenan, 1994; Voris, 2000). The northernmost records of P. textilis in the Northern Territory are also 
indicated. 


Mean between-clade p-distances (± 1 SE) were 0.0236 ± 0.0046 between clade 1 and 2, 0.0290 ± 0.0055 
between clade 1 and 3, and 0.0190 ± 0.0045 between clade 2 and 3. The Milne Bay and Oro haplotypes dif¬ 
fered from each other by a single nucleotide substitution, and from other haplotypes in clade 2 by 4-8 substi¬ 
tutions (p-distance 0.0065-0.0104), and the Merauke haplotype from other clade 3 haplotypes by 1-3 
substitutions (p-distance 0.0013-0.0039). 

Pairwise Wilcoxon signed-ranks tests largely rejected the alternative tree topologies of monophyly of 
either the New Guinea or the Australian populations of P. textilis, although a minority of pairwise compari- 
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sons were insignificant (Table 1). The Shimodaira-Hasegawa tests convincingly rejected the alternative topol¬ 
ogies, and none of the 60004 post-bum-in Bayesian trees was consistent with either of the two alternative 
hypotheses. 


TABLE 1 . Tests of alternative phylogenetic scenarios. 

Wilcoxon signed ranks test Shimodaira-Haseg- Bayesian inference 

awa test 

d(steps) -Z P d(lnL) P proportion of Bayesian trees con¬ 

sistent with alternative hypothesis 

Monophyly of New 10 1.6389-3.1623 0.0016-0.1012 38.90171 <0.001 0/60004 

Guinea haplotypes 

Monophyly of Australian 10 1.8898-2.8868 0.0039-0.0588 42.57724 <0.001 0/60004 

haplotypes 


Discussion 

Biogeography 

The phylogeographic pattern we have identified in this study clearly indicates that the populations of P 
textilis in New Guinea are the result of two separate colonisation events. The populations from the Merauke 
area of Indonesian Papua group with populations of P. textilis from central and central northern Australia, 
whereas the eastern populations from Milne Bay and Oro group with populations from northern Queensland. 
The levels of sequence divergence between the New Guinean haplotypes and their Australian sister groups are 
low in both cases (< 1% p-distance). This suggests that a divergence between the New Guinea populations and 
their Australian sister groups is recent, most likely dating back to the Pleistocene, given most current esti¬ 
mates rates of sequence evolution in squamates (Zamudio & Greene, 1997; Wiister et al., 2002). 

Although there has been much speculation that the New Guinea populations of P. textilis were introduced 
by humans (Slater, 1968; discussed in O’Shea, 1996, Kuch & Yuwono, 2002, and Williams & Wiister, 2005), 
our results are entirely congruent with natural colonisation. We do not consider the abundance or broad distri¬ 
bution of a species to be any indicator of its native or introduced status (Kuch & Yuwono, 2002). In fact, suc¬ 
cessfully introduced species are often found in great abundance across a broad geographic range (e.g. rabbits 
in Australia or the cane toad Bufo marinus in Australia and New Guinea). Indeed, the combination of wide but 
scattered distribution, coupled with apparent rarity, of P. textilis in Central Province and the NCD are more 
suggestive of native status than its abundance elsewhere. Instead, we think strong evidence of the native status 
of P. textilis would be genetic data showing haplotype diversity and divergence: the fact that the haplotypes 
from Oro and Milne Bay Provinces differ by one base pair substitution thus also supports the hypothesis of 
native status of the PNG populations. 

Our data show that the two New Guinea populations possess haplotypes that are nested within the geo¬ 
graphically closest Australian haplotype clades, a pattern that would not necessarily be expected from a 
human introduction. Instead, the pattern observed in P. textilis is consistent with a hypothesis of dispersal dur¬ 
ing times of lowered sea levels during the Pleistocene, which resulted in repeated land connections between 
Australia and New Guinea (Voris, 2000). At various times through the Pleistocene, much of the Sahul Shelf 
was above sea level, leading to more or less wide land bridges between Australia and New Guinea. Voris 
(2000) estimated that a 10 m lowering of sea levels would be sufficient to lead to a narrow land connection 
between Cape York (northern Queensland) and the central southern coast of New Guinea, whereas a lowering 
of 75 m or more would result in extensive emergence of the Sahul Shelf, exposing most of the current sea 
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floor from between Arnhem Land and Aru Island in the west to Cape York and the Gulf of Papua in the east. 
However, during at least part of the Pleistocene, the centre of this shelf was occupied by a large lake, Lake 
Carpentaria, situated within the present-day Gulf of Carpentaria (Torgersen et al., 1983; Keenan, 1994; Voris, 
2000). The presence of this lake essentially meant that there were two dispersal corridors between Australia 
and New Guinea, one along the western side of the Sahul Shelf, from Arnhem Land to southern New Guinea, 
and the other along the eastern side from present-day northern Queensland to the present-day Trans-Fly area 
or the Gulf of Papua, part of which would also have been exposed due to lowered sea levels. Voris (2000) esti¬ 
mated that sea levels dropped to 75m or more below their present level on at least two separate occasions, for 
a total duration of approximately 42,000 years, during the last 250,000 years. Extensive global climatic fluctu¬ 
ations with high-latitude glaciation started occurring at least 2.4 Mya (Shackleton et al. 1984; Hooghiemstra 
& Cleef, 1995), and it therefore seems likely that similar episodes of lowered sea levels have occurred 
throughout the late Pliocene and the Pleistocene, although their timing and extent remain poorly known, and 
are likely to be confounded by tectonic changes in this region. 

Recent genetic studies of several snake taxa indicate that both corridors may have been important in shap¬ 
ing the snake faunas of New Guinea and northern Australia: the Acanthophis rugosus complex clearly used 
the western route, with closely related haplotypes being found in the Merauke area and in the Top End 
(Wiister et al., 2005). On the other hand, water pythons ( Liasis fuscus-mackloti complex—Rawlings et al., 
2004), green tree pythons ( Morelia viridis —Rawlings & Donnellan, 2003) and scrub pythons ( Morelia ame- 
thistina —Harvey et al., 2000) appear to have dispersed solely along the eastern route, although the Liasis fus¬ 
cus-mackloti complex also appears to have undergone separate overwater dispersal from the Northern 
Territory to the lesser Sunda Islands (Rawlings et al., 2004). 

Our analysis of P. textilis reveals the first apparent instance of a species that used both dispersal routes to 
colonise New Guinea from an Australian origin. The close phylogenetic relationship between the eastern PNG 
populations and those from northern Queensland (clade 2) very strongly indicates dispersal along the eastern 
route. Torgersen et al. (1983) suggested that the Fly River drained into Lake Carpentaria rather than the Gulf 
of Papua until the uplift of the Oriomo Plateau during the Pleistocene. As a result, it seems likely that the 
present-day swamplands of the Gulf of Papua coast, where P. textilis and other savanna species such as 
Oxyuranus scutellatus and Pseudechis papuanus do not occur, would have been much drier than is now the 
case. If this hydrographic pattern was still prevailing at the time of the eastern dispersal of P. textilis, then it 
may have facilitated dispersal along the southeastern coast of New Guinea to the eastern tip of the island. 

The Merauke population, on the other hand, is most parsimoniously interpreted as the result of dispersal 
along the western route (Fig. 2). Although P. textilis currently appears to be lacking in the Top End, the north¬ 
ernmost populations from the Northern Territory, presumably assignable to clade 3, are from the base of the 
Peninsula, at Nathan River Station (NT Museum R22316) in the east and at Victoria River Downs (NT 
Museum R3714; see also Gillam, 1979) and “Delamere Road, 2 km S. of Victoria Highway” (15° 20’S, 131° 
37’ E—NT Museum R17151)in the west (Fig. 2). This, plus the fact that the available sampling indicates that 
all the northeastern Queensland populations are assignable to clade 2, suggests that dispersal to the Merauke 
region via the eastern route by clade 3 populations would be an unparsimonious assumption, as it would imply 
repeated exchanges of populations in northern Queensland. More extensive sampling of P. textilis in the 
Northern Territory and northwestern Queensland would provide additional evidence for the most likely dis¬ 
persal routes. 

Within New Guinea, the populations from Central Province and the NCD are situated along the likely dis¬ 
persal route along the southern coast of New Guinea towards the eastern tip of the island. This hypothesis 
therefore predicts that these populations should have haplotypes that group together with those of the Milne 
Bay and Oro specimens. An analysis of additional taxa with wider distributions in northern Australia and 
southern New Guinea (e.g., Morelia spilota, Demansia vestigiata, Furina tristis ) would reveal the relative 
importance of the two dispersal corridors in shaping the snake faunas of Australia and New Guinea. 
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The confirmation of the native status of P. textilis in New Guinea suggests that the paucity of records of 
the species, coupled with a scattered distribution, is most likely a natural phenomenon, and gives little reason 
to predict any rapid change in the relative importance of this species as a snakebite hazard, which would have 
been the case if the species was present as a result of an introduction (Williams & Wiister, 2005). Neverthe¬ 
less, given the adaptability of different species of Pseudonaja in Australia to human-impacted habitats, it is 
possible that the abundance of the species and the frequency of humans being bitten by Pseudonaja in New 
Guinea may increase over coming decades. Indeed, recent fieldwork in New Guinea by the authors suggests 
that both P textilis and Oxyuranus scutellatus may occur in greater densities where savanna has been con¬ 
verted into oil palm monoculture. Similar increased concentrations have also been observed for Micropechis 
ikaheka in New Guinean copra and cocoa plantations (O’Shea, 1994). 

A recent cost-benefit trial of CSL Snake Venom Detection Kits at Port Moresby General Hospital identi¬ 
fied P. textilis as the culprit in four cases from a series of forty envenomed patients (Williams et al., 2006, in 
prep.). In two of these four cases the patients presented 18-24 hours post-snakebite with coagulopathy and 
profound irreversible neurotoxicity and both subsequently died. CSL Brown Snake Antivenom is not pres¬ 
ently registered for clinical use in Papua New Guinea. Continued epidemiological surveillance in New 
Guinea, preferably using immunoassay methods such as the CSL Snake Venom Detection Kits, are therefore 
essential to monitor the situation and ensure that antivenom acquisition and distribution reflects changes in the 
frequency of bites by different species, and thus the needs of bitten patients. Inadequate knowledge and misi- 
dentification of biting snakes has led to medical problems in PNG in the past, when Pseudechis papuanus 
antivenom was stocked and administered to numerous patients, when in reality the vast majority were bitten 
by taipans ( Oxyuranus scutellatus ) (Williams, 2005). 


Taxonomy 

In his description of the first two specimens of P textilis from New Guinea, McDowell (1967) provided basic 
scale and dentition counts for the two specimens, and commented briefly that they “resemble Queensland and 
North Australian material examined, except in having 12, rather than 9 to 11, solid maxillary teeth”. However, 
no indication of the sample size or provenance of the Australian material was provided at the time. Hoser 
(2003) based his description of Pseudonaja textilis pughi primarily on McDowell’s brief description, desig¬ 
nating the specimens described by McDowell as his types, and used the maxillary tooth counts as his main 
diagnostic character. In addition, he briefly described the coloration of living specimens examined by him per¬ 
sonally (without indicating sample size or origin), and explicitly stated that all New Guinea Pseudonaja are 
assignable to his new subspecies. 

Our data reveal that the situation of the New Guinean Pseudonaja is more complex than suggested by 
Hoser (2003). First, the New Guinean P. textilis are diphyletic, with populations from Merauke being 
descended from a different set of Australian populations than the Milne Bay/Oro populations from eastern 
Papua New Guinea. Second, it is also clear that the two New Guinea populations of P. textilis are recent colo¬ 
nisers of New Guinea and share a recent common ancestry with their respective Australian sister populations. 

The extent to which the two New Guinea populations are diagnosably distinct from each other and their 
Australian sister populations remains unclear: McDowell’s (1967) sample size for tooth counts for New 
Guinea was only two, and not stated for Australia (He later stated that his Queensland sample size was “no 
more than eight” [McDowell, pers. comm. 2007 to MOS]). Greer (1997) provided data for an additional 14 
specimens (origin not stated), with a range of 9-11 solid maxillary teeth. However, two of the new specimens 
acquired by the authors (LSUMZ 90636-37) were examined and found to have only 11 solid maxillary teeth, 
demonstrating that the character is not diagnostic for all New Guinea specimens. 
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Hoser’s comments on coloration are also based on an unknown sample size of unknown geographic ori¬ 
gin, and are thus difficult to interpret in view of the documented variability in the coloration and pattern of P. 
textilis from Australia (see Sutherland, 1981, for an example). 

In Table 2, we provide scale counts for most of the known specimens so far documented from Papua New 
Guinea, and for three specimens from Merauke. The scale counts are well within published limits for P. texti¬ 
lis as a whole (Gillam, 1979; Annable, 1985; Cogger, 2000), and there are no clear differences between speci¬ 
mens from eastern PNG and the Merauke area, as far as can be told from our very limited sample. However 
based on the material in Table 2 as well as photographs of four additional Merauke specimens, there appear to 
be consistent differences in colour between the P textilis populations from eastern PNG and from Merauke: 
the Merauke specimens imported to Europe through the Indonesian pet trade are olive or tan to warm mid¬ 
brown, whereas all the PNG specimens examined (including the preserved material in the NMPNG) are dark 
grey-brown to almost black (Fig. 3). 



FIGURE 3. Specimens of the common brownsnake, Pseudonaja textilis, from New Guinea, illustrating colour variation. 
A. AVRU-UPNG PT004, from Heropa plantation, Doboduru, Oro Province, Papua New Guinea. B. AVRU-UPNG 
PT001) from Wamawamana, Milne Bay Province, Papua New Guinea. C. Captive specimen from Merauke, Papua, Indo¬ 
nesia (Klaus Romer, personal collection; photograph by Klaus Romer). Note the warmer colour tone and lighter head 
compared to the PNG specimens. D. NMPNG R 24583, from Port Moresby, NCD, Papua New Guinea. 
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TABLE 2. Scale counts and measurements for known specimens of Pseudonaja textilis from New Guinea. 


Specimen/locality 

Sex 

Dorsals 
(neck, mid¬ 
body, vent) 

Ven- 

trals 

Suhcaudals 

SVL 

TL 

Colour 

LSUMZ 90636: 
Wamawamana, Milne Bay 
Province, PNG 

M 

18-17-15 

205 

>25 divided 

939 mm 

n/a 

(tail incomplete) dark grey-brown 
dorsally with whitish chin and dark 
grey ventral surfaces 

LSUMZ 90637: 
Wamawamana, Milne Bay 
Province, PNG 

M 

18-17-15 

n/a 

62 divided 

764 mm 

152 

(damage to portion of ventral region) 
dark grey-brown dorsally with whitish 
chin and dark grey ventral surfaces 

AVRU-UPNG PT001: 
Wamawamana, Milne Bay 
Province, PNG 

M 

18-17-15 

202 

60 divided 

954 mm 

173 mm 

dark grey-brown dorsally with whitish 
chin and dark grey ventral surfaces 

AVRU-UPNG PT002: 

Heropa, Oro Province, PNG 

M 

17-17-15 

199 

62, divided 

968 mm 

191 mm 

Chocolate brown dorsally, whitish 
chin, grey ventrally 

AVRU-UPNG PT003: 

Heropa, Oro Province, PNG 

M 

17-17-15 

198 

2 single, > 
43 divided 

981 mm 

n/a 

(tail incomplete) Chocolate brown 
dorsally, whitish chin, grey ventrally 

AVRU-UPNG PT004: 

Heropa, Oro Province, PNG 

F 

17-17-15 

201 

63, all 
divided 

977 mm 

180 mm 

Uniformly chocolate brown dorsally, 
whitish chin, grey ventrally 

AMNH 73949 Menapi, Cape 
Vogel, Milne Bay Province, 
PNG 

F 

17-17-15 

-203 

>45 

divided 

-765 

mm 

n/a 

(specimen damaged, tail incomplete, 
no colouration data) 

AMNH 73959: Baiawa, 

Milne Bay Province, PNG 

M 

17-18-15 

205 

>25 

959 mm 

n/a 

(tail incomplete, no colouration data) 

NMPNG 24583: 

Port Moresby, National Capi¬ 
tal District, PNG 

M 

17-17-15 

207 

56, all 
divided 

1195 

mm 

195 mm 

Dark greyish-brown dorsals, chin 
whitish, ventral surfaces grey 

NMPNG 22295: 

Dogura, Milne Bay Province, 
PNG 

M 

17-17-15 

200 

61, all 
divided 

1240 

mm 

212 mm 

Significantly faded 

WW Ptl: 

Merauke, Papua, Indonesia 
(shed skin, courtesy I. Hen- 
drikx) 

M 

17-17-13 

201 

>58, 1-6 
undivided 

n/a 

n/a 

Warm medium brown 

WW Pt2: 

Merauke, Papua, Indonesia 
(shed skin, courtesy Klaus 
Romer) 

M 

17-17-14 

204 

62, 3-4 
undivided 

n/a 

n/a 

Warm medium brown 

Freek Vonk (personal collec¬ 
tion): 

Merauke, Papua, Indonesia 

F 

17-17-15 

206 

63 

960 

193 

Medium brown 


The taxonomic treatment of slightly differentiated allopatric populations has been contentious for many 
years, with different workers using different criteria for the recognition of such isolates as separate species. 
Some have argued that any diagnosably differentiated isolated populations should be regarded as separate spe¬ 
cies, even if in terms of phylogeography they are nested deeply within a widespread species and phylogeo- 
graphically complex species, irrespective of the nature of the difference (e.g., Grismer, 1999). Others have 
retained minimally differentiated recent island populations as part of a single species if phylogeographic pat¬ 
terns place them within more widespread mainland clades (e.g., Pook et al., 2000; Keogh et al., 2005; Castoe 
et al., 2007). Although paraphyly of the mainland populations in gene trees does not preclude recognition of 
the island population as a separate species (Funk & Omland, 2003; Grazziottin et al., 2006), it should stimu¬ 
late a synthetic approach based on a thorough assessment of the systematics of the entire complex. We there¬ 
fore suggest that an assessment of the taxonomic status of the New Guinean populations of P. textilis should 
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involve a comprehensive revision of the entire species complex in both Australia and New Guinea, based, for 
instance, on morphological methods such as multivariate analysis (e.g., Wiister & Thorpe, 1992) in conjunc¬ 
tion with molecular data (e.g., Puorto et al., 2001), rather than a piecemeal approach based on small samples 
and few characters. 
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